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Abstract 



We propose a framework to study optimal trading policies in a one-tick pro-rata 
limit order book, as typically arises in short-term interest rate futures contracts. The 
high-frequency trader has the choice to trade via market orders or limit orders, which 
are represented respectively by impulse controls and regular controls. We model and 
discuss the consequences of the two main features of this particular microstructure: 
first, the limit orders sent by the high frequency trader are only partially executed, and 
therefore she has no control on the executed quantity. For this purpose, cumulative 
executed volumes are modelled by compound Poisson processes. Second, the high 
frequency trader faces the overtrading risk, which is the risk of brutal variations in 
her inventory. The consequences of this risk are investigated in the context of optimal 
liquidation. 

The optimal trading problem is studied by stochastic control and dynamic progra- 
mming methods, which lead to a characterization of the value function in terms of an 
integro quasi-variational inequality. We then provide the associated numerical reso- 
lution procedure, and convergence of this computational scheme is proved. Next, we 
examine several situations where we can on one hand simplify the numerical procedure 
by reducing the number of state variables, and on the other hand focus on specific cases 
of practical interest. We examine both a market making problem and a best execution 
problem in the case where the mid-price process is a martingale. We also detail a high 
frequency trading strategy in the case where a (predictive) directional information on 
the mid-price is available. Each of the resulting strategies are illustrated by numerical 
tests. 



Keywords: Market making, limit order book, pro-rata microstructure, inventory risk, 
marked point process, stochastic control. 



1 Introduction 



In most of modern public security markets, the price formation process, or price discovery, 
results from competition between several market agents that take part in a public auction. 
In particular, day trading sessions, which are also called continuous trading phases, consist 
of continuous double auctions. In these situations, liquidity providers^ continuously set bid 
and ask prices for the considered security, and the marketplace publicly displays a (possibly 
partial) information about these bid and ask prices, along with transactions prices. The 
action of continuously providing bid and ask quotes during day trading sessions is called 
market making, and this role was tradionnally performed by specialist firms. However, due 
to the recent increased availability of electronic trading technologies, as well as regulatory 
changes, a large range of investors are now able to implement such market making strategies. 
These strategies are part of the broader category of high frequency trading (HFT) strategies, 
which are characterized by the fact that they facilitate a larger number of orders being sent 
to the market per unit of time. HFT takes place in the continuous trading phase, and 
therefore in the double continuous auction context, but actual mechanisms that implement 
this general continuous double auction set-up directly influence the price formation process 
and, as a consequence, HFT strategies. 

In this work, we shall focus on the case where the continuous double auction is imple- 
mented by a limit order book (LOB), operated under the pro-rata microstructur^J see |10| 
and pQ. This microstructure can be encountered on some derivatives markets, and espe- 
cially in short-term interest rate (STIR) futures markets, also known as financial futures, 
traded e.g. on LIFFE (London International Financial Futures and options Exchange) or 
on CME (Chicago Mercantile Exchange). We will describe this microstructure in depth 
in Section 2, but the general mechanism of this microstructure is as follows: an incoming 
market order is dispatched on all active limit orders at the best price, with each limit order 
contributing to execution in proportion to its volume. In particular, we will discuss the 
two main consequences of this microstructure on HFT strategies which are the oversizing 
of the best priced slices of the LOB and the overtrading risk. 

Our main goal is to construct an HFT strategy, by means of optimal stochastic control, 
that targets the pro-rata microstructure. We allow both limit orders and market orders in 
this HFT strategy, modelled respectively as continuous and impulse controls, due to consid- 
erations about direct trading costs. From a modelling point of view, the key novelty is that 
we take into account partial execution for limit orders, which is crucial in the pro-rata case. 
For this purpose we introduce a Poissonian model for trades processes, that can be fitted 
to a large class of real-world execution processes, since we make few assumptions about the 
distributions of execution volumes. From a practical trading point of view, we allow the 
HFT to input predictive information about price evolution into the strategy, so that our 
algorithm can be seen as an information-driven HFT strategy (this situation is sometimes 
called HFT with superior information, see [3]). We derive the dynamic programming equa- 
tion corresponding to this mixed impulse/regular control problem. Moreover, we are able to 

In this paper, we call liquidity provider any investor that currently trades with limit orders 
2 This microstructure differs from the price/time microstructure that can be encountered on most cash 
equity markets, see [5]. 
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reduce the number of relevant state variables to one in two situations of practical interest: 
first, in the simple case where the mid-price is a martingale, and second, in the case where 
the mid-price is a Levy process, in particular when the HFT has predictive information 
on price trend, in line with recent studies [5]. We provide a computational algorithm for 
the resolution of the dynamic programming equation, and prove the convergence of this 
scheme. We illustrate numerically the behavior of the strategy and perform a simulated 
data benchmarked backtest. 

High-frequency trading has recently received sustained academic interest. The refer- 
ence work for inventory-based high frequency trading is Avellaneda and Stoikov (2008) [2]. 
The authors present the HFT problem as an inventory management problem and define 
inventory risk as the risk of holding a non-zero position in a risky asset. They also pro- 
vide a closed-form approximate solution in a stylized market model where the controls are 
continuous. Several works are available that describe optimal strategies for HFT on cash 
equities or foreign exchange, e.g. [TT], [7J, [H] or [T2]. Gueant, Tapia and Lehalle ([7]) 
provide extensive analytical treatment of the Avellaneda and Stoikov model. Veraart (|12j) 
includes market orders (that are modelled as impulse controls) as well as limit order in the 
context of FX trading. Guilbaud and Pham ([9]) study market/limit orders HFT strategies 
on stocks with a focus on the price/time microstructure and the bid/ask spread modelling. 
More recently, Cartea, Jaimungal and Ricci ([3]) consider a HFT strategy that takes into 
account influence of trades on the LOB, and give the HFT superior information about the 
security price evolution. A growing literature is dedicated to modelling the dynamics of 
the limit order book itself, and its consequences for the price formation process. A popular 
approach is the Poisson Limit Order Book model as in Cont and de Larrard ([S]). These 
authors are able to retrieve a predictive information on price behavior (together with other 
LOB features) based on the current state of the order book. Finally, in empirical litera- 
ture, much work is available for cash equities e.g. [8], but very few is dedicated to markets 
operating under pro-rata microstructure. We would like to mention the work by Field and 
Large ([6]), which provides a detailed empirical description of pro-rata microstructure. 

This paper is organized as follows: in Section 2, we detail the market model and explain 
the high frequency trading strategy. In Section 3, we formulate the control problem, derive 
the corresponding dynamic programming equation (DPE) for the value function, and state 
some bounds and symmetry properties. We also simplify the DPE in two cases of practical 
interest, namely the case where the price is a martingale, and the case where the investor 
has predictive information on price trend available. In Section 4, we provide the numeri- 
cal algorithm to solve the DPE, and we study the convergence of the numerical scheme 
towards the exact solution, by proving the monotonicity, stability and consistency for this 
scheme. We also provide numerical tests including computations of the optimal policies 
and performance analysis on a simulated data backtest. Finally, in Section 5, we show how 
to extend our model in the optimal liquidation case, i.e. when the investor's objective is to 
minimize the trading costs for unwinding her portfolio. 
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2 Market model 



Let us fix a probability space (O, P, P) equipped with a filtration F = (Pt)t>o satisfying 
the usual conditions. It is assumed that all random variables and stochastic processes are 
defined on the stochastic basis (O, P, F, P). 

Prices in a one-tick microstructure. We denote by P the midprice, defined as a Markov 
process with generator V valued in P. We shall assume that P is a special semimartingale 
such that its predictable finite variation term A satisfies the canonical structure: dA% <C 
d < P >t, with a bounded density process: 

dA t 

(2.1) 



d<P> t 

and the sharp bracket process < P > is absolutely continuous with respect to the Lebesgue 
measure: 

d<P> t = g(P t )dt, (2.2) 

for some positive continuous function g on P. We denote by 5 > the tick size, and we 
shall assume that the spread is constantly equal to 5, i.e. the best ask (resp. bid) price 
is P a := P + | (resp. P b := P — |). This assumption corresponds to the case of the 
so-called one-tick microstructure [6], which can be encountered e.g. on short term interest 
rates futures contracts. 

Trading strategies. For most of investors, the brokerage costs are paid when a transaction 
occurs, but new limit order submission, update or cancel are free of charge. Therefore, 
the investor can submit or update her quotes at any time, with no costs associated to 
this operation: it is then natural to model the limit order strategy (make strategy) as 
a continuous time predictable control process. On the contrary, market orders lead to 
immediate execution, and are costly, so that continuous submission of market orders would 
lead to bankruptcy. Therefore, we choose to model the market order strategy (take strategy) 
as impulse controls. More precisely, we model trading strategies by a pair a = (a make , a take ) 
in the form: 

„make /ro rb\ take ( _ t~ \ 

The predictable processes L a and L b , valued in {0, 1} represent the possible make regimes: 
when LI = 1 (resp. L\ = 1) this means that the investor has active limit orders at the 
best ask price (resp. best bid price) at time t, else, if L" = (resp. L\ = 0) this means 
that the investor has no active order at the best ask price (resp. best bid price) at time t. 
Practical implementation of such rule would be, for example, to send a limit order with a 
fixed quantity, when the corresponding control is 1, and cancel it when it turns to 0. On 
the other hand, (r„) n6 N is an increasing sequence of stopping times, representing the times 
when the investor chooses to trade at market, and £ n , n > are T Tn -measurable random 
variables valued in R, representing the quantity purchased if > or sold if £ n < 0. 

Execution processes in a pro-rata microstructure. The pro-rata microstructure (see 
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|10| for extensive presentation and discussion) can be schematically described as follows^] 
when a market order comes in the pro-rata limit order book, its volume is dispatched 
among all active limit orders at best prices, proportionnally to each limit orders volumes, 
and therefore create several transactions (see Figure IT]). 




Figure 1: Schematic view of the pro-rata market microstructure. 

This pro-rata microstructure fundamentally differs from price-time microstructure |9j for 
two reasons: first, several limit orders at the best prices receive incoming market order flow, 
regardless of the time priority, and second, market makers tend to oversize their liquidity 
offering (that is, posting limit order with much higher volume than they actually want to 
trade) in order to increase their transaction volume. For example, on the three-months 
EURIBOR futures contracts, the liquidity available at the best prices is 200 times higher 
than the average transaction size. Therefore, to model a market making strategy, one must 
take into account the fact that limit orders are always oversized, so that the executed volume 
is a random variable on which the market maker has no controQ More precisely, let N a 
(resp. N b ) be a Poisson process of intensity A a > (resp. A 6 ), whose jump times represent 
the times when execution by a market order flow occurs at best ask (resp. best bid), and we 
assume that N a and N b are independent. Let (Cn)«.eN* and (Cn)™eN* be two independent 
sequences of i.i.d integrable random variables valued in (0, oo), of distribution laws [i a and 
fjr, which represent the transacted volume of the n th execution at best ask and best bid. 
We denote by v a {dt,dz) (resp. v b (dt,dz)) the Poisson random measure associated to the 
marked point process (-/V a , (Cn)neN*) (resp. (N b , (Cn)neN*)) °f intensity measure \ a fi a (dz)dt 

3 For a detailled description of actual trading rules, and a general overview of STIR futures trading, we 
refer to [T] and references therein. 

4 This differs from the price-time microstructure case, in which the market maker can control the upper 
bound of the execution volume by adjusting the limit order volume. 
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(resp. X b fi b (dz)dt), which is often identified with the compound Poisson processes 



N? 



0? 



n=l -° J0 n=l 

representing the cumulative volume of transaction at ask, and bid. Notice that these 
processes model only the trades in which the investor has participated. 

Cash holdings and inventory. The cash holdings process X and the cumulated number 
of stocks Y (also called inventory) hold by the investor evolve according to the following 
dynamics: 



JV<> 



= / / Z^(dt,dz), 4 = J2< b n= I I ZU b (dt,dz). (2.3) 



JO 



dX t 
dY t 

Y —Y 

1 T„, 1 T„ — 



L a t {p t - + -)dr t -L b t {p t 
L a t dr t , 



-)d#\, T n <t< T n+1 



L° t d®l 



T n <t < T n+ i 



,5 



<nPr n - l&Kg +e) -e , 



in 



(2.4) 
(2.5) 
(2.6) 
(2.7) 



The equations (2.4)-(2.5) model the evolution of the cash holdings and inventory under 



a limit order (make) strategy, while equations (2.6)- (2.7) describe the jump on the cash 



holdings and inventory when posting a market order (take) strategy, subject to a per share 
fee e > and a fixed fee £q > 0. In the sequel, we impose the natural admissibility condition 
that the size of the market order should not be larger than the current inventory, i.e. |£ n | < 
\Y Tn — |, n > 0, and we shall denote by A the set of all admissible make and take strategies 
a = (a make ,a take ). 



Remark 2.1 Let us define the process Vt = Xt + YfPt, which represents at time t the 



marked-to-market value of the portfolio (or book value of the portfolio). From (2.4)-(2.5)- 



(2.6)-(2.7), we see that its dynamics is governed by: 

5 



dV t 



(L b dO b + L?d^) + Y t -dP t , 



-I£«l(2+ e ) 



£ . 



(2.8) 
(2.9) 



In equation (2.9), we notice that a trade at market will always diminish the marked to 



market value of our portfolio, due to the fact that we have to "cross the spread" , hence trade 



at a least favorable price. On the other hand, in equation (2.8), the term J ^(L b d , d b -\-Lfd , &f) 
is always positive, and represents the profit obtained from a limit order execution, while 
the term fY t -dPt represents the portfolio value when holding shares in the stock, hence 
inducing an inventory risk, which one wants to reduce its variance. 



3 Market making optimization procedure 
3.1 Control problem formulation 

The market model in the previous section is fully determined by the state variables (X, Y, P) 
controlled by the limit/market orders strategies a = (a make , a take ) 6 ^4. The market maker 
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wants to optimize her profit over a finite time horizon T (typically short term) , while keeping 
control of her inventory risk, and to get rid of any risky asset by time T. We choose a mean- 
variance optimization criterion, and the goal is to 



maximize E 



X T -1 Y t ±d<P> t 



over all strategies a £ A, s.t Yp = 0. (3.1) 



The integral J Y f 2 _d < P > t is a quadratic penalization term for holding a non zero 
inventory in the stock, and 7 > is a risk aversion parameter chosen by the investor. 

Let us now rewrite problem (3.1 ) in a more standard formulation. Notice indeed that one 
can remove mathematically the constraint Yt = on the inventory control, by introducing 
the liquidation function: 
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x + yp- \y\(- +e) - e , 



which represents the cash obtained after an immediate liquidation of the inventory via a 
market order. Thus, problem (3.1) is formulated equivalently as 

cT 



maximize E 



L(X T , Y T , P T )-1 I Y t 2 g{P t )dt over all strategies a £ A, (3.2) 







sup E tjX ^ p L(X T , Y T , P t ) - 7 / Yf e (P s )ds 



(3.3) 



Here, given a £ A, Et tX ,y,p denotes the expectation 



where we used also (2.2). Let us then define the value function for the problem (3.2): 

v{t,x,y,p) = 

for t G [0,T], (x,y,p) G M 2 
operator under which the process (X, Y, P) solution to (2.4)-(2.5)-(2.6)-(2.7) with initial 
state (X t - , Y t - , P t - ) = (x,y,p), is taken. Problem (3.3) is a mixed impulse/regular con- 
trol problem in Markov model with jumps that we shall study by dynamic programming 
methods. 

First, we state some bounds on the value function. 

Proposition 3.1 There exists some constant Kp (depending only on the price process and 
j) such that for all (t,x,y,p) £ [0,T] xtfxF, 



L(x,y,p) < v(t,x,y,p) < x + yp + -(\ a fl a + X b ff)(T - t) + K P , (3.4) 
where j2 a = J °° z\i a {dz), fl b = z\${dz) are the mean of the distribution laws \i a and fJ 3 . 



Proof. The lower bound in (3.4) is derived easily by considering the particular strategy, 
which consists of liquidating immediately all the current inventory via a market order, and 
then doing nothing else until the final horizon. Let us now focus on the upper bound. 
Observe that in the definition of the value function in (3.3), we can restrict obviously to 
controls a £ A s.t. 



E 



Y t z d <P> t 



< 00. 



(3.5) 
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For such strategies, we have: 



E 



t,x,y,p 



< E 



t,x,y,p 



L{X T ,Y T ,P T )- 1 J Y s 2 d<P> s 
V T -1 Y s 2 d < P > s 



< x + yp + E tiX ^ p (r T _ t + + ^ Y s - dP s - 7 ^ Y 2 d < P > £ 

= x + y P + V t \ 5 (r T _ t + $ b T _ t )+ f T (Y s -9 s - 7 Y s 2 )d<P>. 



Here, the second inequality follows from the relation ( 2.8 ), together with the fact that L a ,L 
< 1, i9 a , t? 6 are increasing processes, and also that jumps of V are negative by (2.9). The 



last equality holds true by (2.1 ) and the fact that J Y_dM is a square- integrable martingale 



from (3.5), where M is the martingale part of the semimartingale P. Since 6 is bounded 



and 7 > 0, this shows that for all strategies a satisfying (3.5), we have: 

rT 



E ttX , y>p \L(X T , Y T , P T ) - 7 J Yfd < P > s 

< x + yp + 6 -E[r T _ t + $ b T _ t } + KE[<P> T }, 

for some positive constant K, which proves the required result by recalling the character- 
istics of the compound Poisson processes # a and i? 6 , and since < P >t is assumed to be 
square-integrable. □ 



Remark 3.1 The terms of the upper bound in (3.4) has a financial interpretation. The 



term x + yp represents the marked-to-market value of the portfolio evaluated at mid-price, 
whereas the term Kp stands for a bound on profit for any directional frictionless strategy 
on the fictive asset that is priced P. The term T l (\ a ji a + \ b Jl h ){T — t), always positive, 
represents the upper bound on profit due to market making, i.e. the profit of the strategy 
participating in every trade, but with no costs associated to managing its inventory. 



3.2 Dynamic programming equation 

For any (£ a ,£ b ) G {0, l} 2 , we introduce the non-local operator associated with the limit 
order control: 



where 



T a (f)(t,x,y,p) = X a [<f)(t,x + z(p + -),y- z,p) - <P(t,x,y,p)]n a (dz) 
Jo z 

f°° 5 
T b (j)(t,x,y,p) = X b [<f)(t,x- z(p- -),y + z,p) -<p(t,x,y,p)]fi b (dz), 



(3.6) 
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for (t, x, y,p) [0,T] xlxRxP. Let us also consider the impulse operator associated with 
admissible market order controls, and defined by: 



M<fi(t,x,y,p) 



sup 

ee[-\y\,\y\\ 



(t,x-ep- |e|(- +e) -e ,y + e,p). 



The dynamic programming equation (DPE) associated to the control problem (3.3) is 
a quasi-variational inequality (QVI) in the form: 



mm 



Ov 
di 



~pa no 

sup L ' v + 75 , v — M.v 

(£ a ,e b )e{o,i} 2 



together with the terminal condition: 

v(T, .) 



L, on 



0, on [0,T) x R 2 x P, (3.7) 



(3.8) 



where we denoted by g the function: g(y,p) = y 2 g(p)- This DPE may be written explicitly 
as: 

dv 



mm 



ot 



Vv 



/ f°° 5 \ 

A a (J [v(t,x + z(p+-),y-z,p) -v(t,x,y,p)]fi a {dz)) + (3.9) 



- A 6 ( [v(t,x-z(p- ~),y + z,p) - v(t,x,y,p)]fi b (dz)J + + jy 2 g(p) 

v(t,x,y,p)- sup v(t,x-ep-\e\(~ + e)-e Q ,y + e,p) 
ee[-\y\,\y\] 2 

for (t,x,y,p) G [0,T) x M? x P, together with the terminal condition: 

v(T,x,y,p) = x + yp- \y\(- + e) - e , V(x, y,p) G R 2 x P. 



(3.10) 



By standard methods of dynamic programming, one can show that the value function in 
(3.3) is the unique viscosity solution under growth conditions determined by (3.4) to the 
DPE (5.2)-(3.10) of dimension 3 (in addition to the time variable). 



3.3 Dimension reduction in the Levy case 

We now consider a special case on the mid-price process where the market making control 
problem can be reduced to the resolution of a one-dimensional variational inequality invol- 
ving only the inventory state variable. We shall suppose actually that P is a Levy process 
so that 

"P/p = c p , and g is a constant, (3-H) 

where ip is the identity function on P, i.e. Ip(p) = p, and g > 0, c p are real constants 
depending on the characteristics triplet of P. Two practical examples are: 

• Martingale case: The mid-price process P is a martingale, so that Vfy = 0. This 
martingale assumption in a high-frequency context reflects the idea that the market maker 
has no information on the future direction of the stock price. 
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• Trend information: To remove the martingale assumption, one can introduce some 
knowledge about the price trend. A typical simple example is when P follows an arithmetic 
Brownian motion (Bachelier model). A more relevant example is described by a pure jump 
process P valued in the discrete grid 5Z with tick 5 > 0, and such that 



P{P t+ h-Pt = S\Tt) 
P(P t+h -P t = -6\T t ) 
P(\P t+h -Pt\>S\Tt) 



ir + h + o(h) 
ir~h + o(h) 
o(h), 



where 7r + , 7r~ > 0, and o(h) is the usual notation meaning that lim^o o(h)/h = 0. Relation 

vj5, where w = tt + — ir~ represents a constant information 
(ir + + 7r~)<5 2 . In a high-frequency context, this model is 



(3.11) then holds with c F 



about price direction, and g 
of practical interest as it provides a way to include a (predictive) information about price 
direction. For example, work have been done in [5] to infer the future prices movements 
(at the scale of a few seconds) from the current state of the limit order book in a Poisson 
framework. In this work, as well as in our real data tests, the main quantities of interest are 
the volume offered at the best prices in the limit order book, also known as the imbalance. 



mm 



0, on [0,T) x 



In this Levy context, the value function v is decomposed into the form: 

v(t,x,y,p) = L(x,y,p) + w(t,y), 
where w is solution to the integral variational inequality: 

yc p + jgy 2 — I a w — I b w , w — M.w 
together with the terminal condition: 

w(T,y) = 0, VyeM, 
where X a and X b are the nonlocal integral operators: 

l a w(t,y) = X a 



dw 
~dt 



(3.12) 

(3.13) 
(3.14) 



s s 

w(t,y- z) - w(t,y) + z- + (- + e)(\y\ 



l b w(t,y) 



A 



8 $ 

w{t, y + z)- w(t, y) + z- + (- + e)(\y\ 



and M is the nonlocal operator: 



Mw(t,y) 



sup 

ee[-\y\,\y\] 



w(t,y + e) - (- +e)(\y + e\ + |e| 



y-z\)]fi a (dz)) 
y + z\)]^ h {dz 



\y\) -^o 



The interpretation of the decomposition (3.12) is the following. The term L(x,y,p) rep- 
resents the book value that the investor would obtain by liquidating immediately with a 
market order, and w is an additional correction term taking into account the illiquidity 
effects induced by the bid-ask spread and the fee, as well as the execution risk when sub- 
mitting limit orders. Moreover, in the Levy case, this correction function uu depends only 
on time and inventory. From (|3.4|), we have the following bounds on the function w. 



< w(t,y) <(^+e)\y\ + ^(\ a p, a + \" fl ' , HT-l) + K,>. V(/.//) e [(.)..'/'] x 



Actually, we have a sharper upper bound in the Levy context. 
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Proposition 3.2 Under (3.11), we have: 



< w(t,y) < (T-t) 

for all (t, x, y, p) G [0, T] x I 2 x P. 
Proof. For any (x,y,p) 6 M 2 x 



u ^ + A a (5 + e)/2 a + A b ((5 + e)^ 



we notice that 



L(x,y,p) 
eo + (2 +e) 



sup L(x - ep - \e\{- + e) - e ,y + e,p) 



e£[-\y\,\y 

\y\ + 



inf |e| + |y + e| 
e£{-\y\,\y\] 



E > 0. 



(3.15) 



We also observe that for all z > 0: 



L(x + ^),y~ z,p) - L(x,y,p) 



z^ + h + \y-A) 



2 y 2 
< {S + e)z 



and similarly: 



L(x - z(p - -),y + z,p) - L(x,y,p) < (5 + e)z. 



(3.16) 



(3.17) 



Let us then consider the function (f>(t, x, y,p) = L(x, y,p) + (T — t)u, for some real constant 
u to be determined later. Then, 4>(T, .) = L, and by (3.16 )-( 3.17| ), we easily check that: 



dt 



sup L ' <p + jg 



> u - X a (5 + e)fl a - X b (5 + e)fl b - yc P + iy 2 p. 

The r.h.s. of this last inequality is a second order polynomial in y and therefore it is always 
nonnegative iff: 

c 2 P -A 1P (u- X a (5 + e)fl a - X b (5 + e)fl b ) < 0, 

which is satisfied once the constant u is large enough, namely: 

r 2 

u > u := i^ + X a (5 + e)fl a + X b (5 + e)fl b . 
47/9 

For such choice of u = u, and denoting by <f> the associated function: (p(t,x,y,p) = 
L(x,y,p) + (T — t)u we have 



dt 



sup C 

(e a ,e b )£{o,i} 2 



£ a l b 



■ + 19 > 0, 



which shows, together with (3.15), that eft is a supersolution of (3.7)-(3.8). From comparison 
principle for this variational inequality, we deduce that 

v < 4> on [0, T] x M 2 x P, 

which shows the required upper bound for w = v — L. □ 



Finally, from (3.13)-(3.14), and in the case where A a = A b , [i a = fj, b , and by stressing 
the dependence of i« in c p , we see that w satisfies the symmetry relation: 

w(t,y,c p ) = w(t,-y,-c P ), V(t, y) & [0, T] x R. (3.18) 
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4 Numerical resolution 



In this section, we focus on the numerical resolution of the integral variational inequality 



(3.13)-(3.14), which characterizes the reduced value function of the market-making problem 



in the Levy case. 



4.1 Numerical scheme 



We provide a computational scheme for the integral variational inequality (3.13). We first 



consider a time discretization of the interval [0, T] with time step h = T/N and a regular 

time grid IV = {tk = kh , k = 0, . . . , N}. Next, we discretize and localize the inventory 

M 

state space on a finite regular grid: for any M > and Ny 6 N, and denoting by Ay = — — , 

Ny 

we set: 



{ yi = iAy, i = -Ny,...,Ny}. 



We denote by Proj M (y) := — M V (y A M), and consider the discrete approximating distri- 
bution of /i a and fx b , defined by: 



A a = Yl /A^Ay;(i + l)Ay))^ A , 



A 6 = /A[*Ay;(* + l)Ay))5 iAy , 



with 5 X the Dirac measure at x. We then introduce the operator associated to the explicit 



time-space discretization of the integral variational inequality (3.13): for any real- valued 
function ip on [0, T] x R, i £ [0, T], and y G R, we define: 



S^ M (t f y,<p) = max T h < A r> M (t,y,<p) ; M h ^' M (t,y^) 



where 



T hAY ' M (t,y,<p) = ip(t,y)-h 70 y 2 + hycp 

+ X a hU [p(t,Vm] M (y - z)) - <p(t,y)]fi a (dz) 



and 



+ J [\z + { 5 -+e)(\y\-\y-z\)]^{dz))_ 

roc 

+ X b h( y J [<p(t,PToj M {y + z)) - <p(t,y)]fi\dz) 

[^z+d + E)(\y\-\y + z\)]»\dz))_ 



M 



'(t,y,<p) 

5 

sup [ip(t,Proi M (y + e))-(- + e){\y + e\ + \e\-\y\)-e ]. (4.1) 



eeY M n[-|y|,|y|] 
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By recalling that x + = max( g { 0i i} fa, we see that the operator T n ' LX) 
also as: 



T hA ^ M (t,y,ip) 



-hjgy + hycp + max 

e a ,£ b e{o,i} 



+ \ a hl° 



" A> xI may be written 
(4.2) 



<p(t,y)(l - X a ht - X b h£ b ) 
¥>(*,Proj M (y - z))(i a (dz) 



o 



f°° S 

+ J o [-z + (- + e)(\y\-\y-z\)}^(dz: 

oo 



+ X b hf b U <p{t,Vro] M {y + z))iL h {dz) 



o 

+ jf [^ + (|+^)(|y|-|y + ^|)]A^) 



Note that the integral terms involving <p in S h '^ Y ' M (t,y, (p) are in fact finite sums, and 
therefore are readily computable. We also assume, for simplicity sake, that the terms 
io°° & 2 Z + (I + £ )(lyl - If - ^IK(^) and J °° §« + (§ + e){\y\ - \y + z\)^ b (dz) are exactly 
computable. 



We then approximate the solution w to (3.13)-(3.14) by the function w h '^ Y ' M on Tjv X 
Yj\,/ solution to the computational scheme: 



h,Ay,M / + 



w h,Ay,M 



(tk,y) 



(4.3) 
> M (t k+1 ,y,w h ' AY > M ) , A; = 0,...,iV-l,yGY M . (4.4) 



4.2 Convergence of the numerical scheme 



In this section, we study the convergence of the numerical scheme (4.3)-(4.4) by showing the 
monotonicity, stability and consistency properties of this scheme. We denote by C^([0, T] x 
K) the set of bounded continuously differentiable functions on [0, T] x R, with bounded 
derivatives. 

Proposition 4.1 (Monotonicity) 

For any h > s.t. h < ^ a ^ ^ e operator S h,Ay ' M is non- decreasing in ip, i.e. for any 
(t,y) € [0,T] x R and any ip,ip G C & X ([0,r] x R) , s.t. (p<tp : 



Proof. From the expression (4.2), it is clear that T h ' Y ' M (t, y, ip), and then also S h > Y ' M (t, y, 



is monotone in <p once 1 — X a h — X h > 0. 



□ 



Proposition 4.2 (Stability) 

For any h,Ay,M > there exists a unique solution w h ' AY ' M to (4.3)-(4.4), and the se- 
quence (w h ' Y ' M ) is uniformly bounded: for any (t,y) E T/v x Ym, 



r ^2 



e < ^' Ay ' M (t,y) < (T - t) + A a (<5 + e)fi a + A (<5 + e)p 

L4 7 p 
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Proof. Existence and uniqueness of w h,/ ^ Y,M follows from the explicit backward scheme 



(4.3)-(4.4). Let us now prove the uniform bounds. We consider the function 

(T-t) 



j?- + X a (8 + e)fl a + X\5 + e)fl b 



and notice that H!*(t) > S h > Y,M {t + h,y,^*) by the same arguments as in Proposition 
Moreover, we have, by definition, w h ' AY,M (T,y) = ty*(T) = 0, and therefore, a direct 



3.2 



recurrence from gjk-d44k shows that w hAY ' M (t, y) < **(t) for all (t, y) E T n xY M , which 
is the required upper bound for w h ' AY ' M . 

On the other hand, we notice that S h ' Y,M (t, 0, (p) > <p(t,0) for any function ip on 
[0, T] x R, and t E [0, T], by considering the "diffusive" part of the numerical scheme with 
the particular controls i a = £ b = 0. Therefore, since w hAY ' M (T, 0) = 0, we obtain by 



induction on (4.3)-(4.4) that w ' Y > (t,0) > for any t E Tjv- Finally, considering the 



□ 



obstacle part of the numerical scheme, with the particular control e = —y, shows that 
w h,AY,M(t,y) > w h ' AY ' M (t,0) - e > -e for any (t,y) E IV X Y M , which proves the 
required lower bound for yj hAr,M . 

Proposition 4.3 (Consistency) 

For all (t,y) E [0, T) xl and tp £ C£([0,T] x R), we have 



lim — 

(h, Ay , M) -> (0, 0, oo) /l 
(*',„') -> (t lV ) 



(t, y) - yc P + 16 y 2 - l a p(t, y) - 7*<p(t, y) 



(4.5) 



and 



lim M hAv ' M (l! + h,y',ip) = M<p(t,y) 

(h, A Y .M) -> (0, 0, oo) 
(*',»') -> (t,v) 



(4.6) 



Proof. The consistency relation (4.6) follows from the continuity of the function (t,y,e) 
6 — 

— > tp(t,y + e) — (- + e)(\y + e| + |e| — |y|) — eo- On the other hand, we have for all (t',y') 
E[0,T)xM, 

1 

hi 



ip{t',y')-T^ M {t' + h,y'M 



[<p(t', y') - <p(t' + h, y')] - y'c P + 1P y' 2 (4.7) 



where 



I^ A - M (t,y,^) = A°(y [p(t,Proj M (y-z))-^,y)]£ a (dz) 



+ y [^ + (2+^)(lyl-ly-^l)]^ a (^))_ 

/-OO 

r (t,y,V9) = A ft ^y [^(t,ProjV(y + 2))-y(t,y)]A 6 (<te) 

+ J™ [ 5 -z + ( 5 -+e)(\y\-\y + z\)}v b (dz))_ 
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The three first terms of (4.7) converge trivially to — —-(t,y) — yc p + jgy 2 as h goes to 



zero and (t 1 , y') goes to (t, y). Hence, in order to get the consistency relation, it remains to 
prove the convergence of Ia ,Av,M (t' + h, y' , <p) to I a ip(t, y) as (h, Ay, M) goes to (0, 0, oo), 
and (t' , y') goes to (t, y) (an identical argument holds for 1^' Y,M (t 1 + h, y' , (p)). By writing 
that \x + — x',\ < \x — x'\, we have 



i^Y> M (t' + h,y',<p)-iMt,y) 

< \ a \<ptf + h,y')-cp(t,y)\ 



+ A a 



<p(t' + h,Proj M (y' - z))jl a (dz) 



< \ a \<p(t' + h,y')-<p(t,y)\ 



+ A a 



+ A c 



cM- 



(p(t' + h,y' - z)fi a (dz) 



cM- 



<p(t,y- z)n a (dz) 



<p(t,y- z)ij, a (dz) 



ip(t' + h,-M)fi a (dz) - / <p(t,y - z)fi a (dz) 

lM+y' JM+y' 

< \ a \<p(t' + h,y')-<p(t,y)\ 

/•oo 

+ A a / \cp(t' + h,y' -k(z)) - ip(t,y-z)\fi a (dz) 



+ 2A a |M| 00 ^ a ([M + 2 /,oo)), 
where we denote by n(z) = \_— — J Ay. Here [z\ denotes the largest integer smaller than z. 

Ay 

Since the smooth function ip has bounded derivatives, say bounded by \\(p^ \\oo, it follows 
that 



lZ AY > M (t' + h,y',<p)-l a <p(t,y) < X a \\p {1) \\oo{h + 2\y'-y\+A Y ) 



+ 2A a |M| oo/ / l ([M + 2 /,oo)), 



which proves that 



rl ' Av ' M (t' + M',^) = l a p(t,y), 
hence completing the consistency relation (4.5). 



lim 1^ Y ' 

(h, A Y ,M) -+ (0, 0, oo) 
(t',y') -> (t,y) 



□ 



Theorem 4.1 (Convergence) 

The solution w h > Y,M to the numerical scheme ('(4.3)-(4.4) y ) converges locally uniformly to 
w on [0,T) x R, as (h,A Y ,M) goes to (0,0, oo). 

Proof. Given the above monotonicity, stability and consistency properties, the convergence 
of the sequence (w h > Y > M ) towards w, which is the unique bounded viscosity solution to 



(3.13)-(3.14), follows from [3]. We report the arguments for sake of completeness. From 
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the stability property, the semi-relaxed limits: 

w*(t,y) = 



liminf ic 

(h, Ay, M) -> (0, 0, oo) 
(t'.y') -> (t,») 



lim sup 



w 



h ' AY ' M {t',y'), 
hAY ' M (t>,y>), 



(h, Ay, Af) -» (0, 0, oo) 

(t',j/') -> (t,») 



are finite lower-semicontinuous and upper-semicontinuous functions on [0, T] xl, and inherit 
the boundedness of (w h, ^ Y,M ). We claim that w* are w* are respectively viscosity super and 
subsolution of (3.13)-(3.14). Assuming for the moment that this claim is true, we obtain 
by the strong comparison principle for (3.13)-(3.14) that w* < w*. Since the converse 
inequality is obvious by the very definition of w* and w*, this shows that w* = w* = w 
is the unique bounded continuous viscosity solution to (3.13)-(3.14), hence completing the 
proof of convergence. 

In the sequel, we prove the viscosity supersolution property of w* (a symmetric argument 
for the viscosity subsolution property of w* holds true). Let (t, y) S [0,T)xE and p a test 
function in C^([0,T] x R) s.t. (t, y) is a strict global minimimum point of — tp. Then, 
one can find a sequence (t' n , y' n ) in [0, T) x E, and a sequence (h n , Ay, M n ) such that: 



(t,y), (h n ,A Y ,M n 



(0,0,oo), w 



w* (t,y), 



(t' n , y' n ) is a global minimum point of w 



/in, AS- 



cp. 



Denoting by ( n = (w hn ' A Y ,M " — p)(t' n , y' n ), we have w h nAY> M ™ > ip+( n . From the definition 
of the numerical scheme 5' l " ,A y ,M " , and its monotonicity, we then get: 



( n + <p(t' n ,y' n ) = w 



(^n> Vn) 



> 



S h ^ M -(t' n + h n ,y' nli p + Cn) = S h "^ M "(t' n + h n ,y' n ,p) + C n 



max 



T KA^M n{t , n + K ^ nM > M h ^ M -{t' n + K,y' n M 



h n Av,M n 



which implies 



mm 



V 



tin Ay, M, 



'(t n + h n ,y' n ,ip) 



hr, 



, <p(t' n ,y' n )-M h ^v M "(t' n + h n ,y' n ,<p) 



> 0. 



By the consistency properties (4.5)-(4.6), and by sending n to infinity in the above inequa- 
lity, we obtain the required viscosity supersolution property: 



mm 



dip 
~dt 



(t, y) - yc P + 1Q y 2 - I a p(t, y) - l b p(i, y) , tp(t, y) - M<p(t, y) 



> 0. 



□ 



4.3 Numerical tests 

In this section, we provide numerical results for the (reduced-form) value function and 
optimal policies in the martingale case and the trend information case, and a backtest on 
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Parameter 


Value 


6 


12.5 EUR/contract 


e 


1.05 EUR/contract 


£o 





A 


0.05s- 1 


1* 


exp(l/A) 


P> 


20 contracts 


7 


2.5. 10" 5 


T 


100 s 



Parameter 


Value 


Ny 


100 


N T 


500 


N m 


20 



(b) Discretization parame- 
ters 



(a) Market and risk parameters 

Figure 2: Parameters for numerical results. 



simulated data for the trend information case. Parameters for these numerical tests are 
shown in Figure [2] 

This set of parameters are chosen to be consistent with calibration data on the front 
maturity for 3- months EURIBOR future, see for example [6]. Within this section, and for 
this set of parameters, we will denote by w h the value function and by a* the make/take 



strategy associated with the backward numerical scheme (4.3)-(4.4). Given a generic con- 
trolled process Z and a control a G A, we will denote Z a the process controlled by a. Unless 
specified otherwise, such processes will be supposed to start at zero: typically, we assume 
that the investor starts from zero cash and zero inventory at date t = in the following 
numerical tests. Finally, we will write indifferently w h (t,y, cp) or w h {t,y) := w h (t,y,0) to 
either stress or omit the dependence in cp. 



• The martingale case: in the martingale case, we performed the algorithm (4.3)-(4.4) 
with parameters shown above. 

Figure [3] displayed the reduced-form value function w h on [0, T] x [— Ny ; Ny] ■ This result 



illustrates the linear bound (3.2) as noticed in proposition 3.2, and also the symmetry of 
w h as pointed out in ( |3.18 ) . We also observe the monotonicity over M + and M~ of the value 
function w h (t, .). 

In Figure |4j we display the optimal make and take policies. The optimal take policy 
(on the left side) is represented as the volume to buy or sell with a market order, as a 
function of the time and inventory (t,y) G [0, T] x [— Ny,Ny]. We notice that a market 
order only occurs when the inventory becomes to large, and therefore, the take policy can 
be interpreted as a "stop-loss" constraint, i.e. an emergency rebalancing of the portfolio 
when the inventory risk is too large. 

The optimal make policy is represented as the regime of limit orders posting as a function 
of the time and inventory (t, y) G [0, T] x [—Ny; Ny]. For sake of simplicity, we represented 
the sum of £ a and & on the map. The meaning of this surface is as follows: means that 
there is no active limit orders on either sides, 2 means that there is active limit orders on 
both bid and ask sides, and 1 means that there is only one active limit order either on 
the bid or the ask side, depending on the sign of y (if y < only the bid side is active, 
and if y > only the ask side is active). We notice that when close to maturity T, the 
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KKI 

Figure 3: Reduced form value function w h . 




(a) Optimal take policy. (b) Optimal make policy 



Figure 4: Numerical results for the martingale case: representation of optimal make and 
take policies a*. 

optimal strategy tends to be more agressive, in the sense that it will seek to get rid of any 
positive or negative inventory, to match the terminal liquidation constraint. Moreover, we 
notice that close to date 0, the dependence in t seems to be negligible, which indicates that 
a "stationary regime" may be attained for large T. 

• The trend information case: in this case, we provide a backtest of the optimal strategy 
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on simulated data in addition to the plots of the value function w h and optimal policy a*. 







-50 





Inventory 



50 



100 



-50 



50 
Inventory 



100 



(a) Value function w h at date t — 0. (b) Optimal policy a* at date t = 0. 

Figure 5: Value function and optimal policy for the trend information case. 

Figure [5] displays the value function and optimal policy at date t = 0, in the plane 



(y,cp). The value function has central symmetry properties as expected in (3.18), and 



should be read as follows: dark green zones represent situation where a market order to 
buy must be sent, light green means that a limit order is active only at bid, white means 
that limit orders are active on both sides, light red means that a limit order is active only at 
ask, and dark red means that a market order to sell must be sent. The value function also 
increases with \cp\. This effect can be interpreted as the gain in performance due to the 
superior information on price trend cp. The interpretation of this extra performance due to 
cp is that the optimal policy avoids part of the adverse selection risk by using this predictive 
information about price movements. Let us provide a qualitative example: assume that 
after the high frequency trader acquired a positive inventory, the adverse selection effect 
implies that price should go down; therefore, using the fact that in this case we should have 
cp < 0, the optimal strategy will be either to cancel the bid limit order (light red zone) 
and keep ask limit order active, or depending on the value of \cp\, send a market order to 
get rid of our positive inventory (dark red zone). 

We performed a benchmarked backtest on simulated data and a performance analysis 
in this case. The benchmark strategy is made of constant controls (a.k.a symmetric or 
constant strategy): 

benchmark t _ make , benchmark 



a 

, make , benchmark 



a 



(1,1) 



\0) 



In order to make our simulated data backtest closer to the reality, we chosed to slightly 
deviate from the original price model, and use a varying price trend. We simulate a price 
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process model given by 



P t = P + 6(N+-Nf), 



where N + and iV are the Euler scheme simulation of Cox processes of respective intensities 
7r + and 7r~ defined as follows 



dir. 



7T + 7T 



= K = g/5 2 

= —Owtdt + odB t 



where K > 0, 6 > and a > are positive constants, and B is an independent Brownian 
motion. Note that we choosed the sum ir + + ir~ to be the constant K, for simplicity sake: it 
means that, disregarding the direction of price variation, the mean number of price change 
per second is assumed to be constant P {\Pt+h — Pt\ = $) = Kh + o(h), which provides an 
easy way to calibrate the parameter K while reducing the dimension of the simulation. The 
interpretation of this simulation model is as follows: we add an exogenous risk factor B, 
which drives the price trend information w as an Ornstein-Uhlenbeck process. Notice that 
this supplementary risk factor B is not taken into account in our optimization procedure 
and thus has a penalizing impact on the strategy's performance: therefore it does not spoil 
the backtest. This model choice for the process (zut) is an convenient way to simulate the 
real-world situation, where the high-frequency trader continuously updates her predictive 
information about short-term price movements, based e.g. on the current state of the limit 
order book. Therefore, qualitatively speaking, our optimization procedure is consistent with 
this simulation model if we choose 6 and a s.t. the variation of the (reduced-form) value 
function w due to predictive information is very small compared to the variation of the value 
function due to other market events (e.g. an execution event). This assumption is consistent 
with HFT practice since the HF trader is able to adapt very quickly to a modification of 
this predictive information. Backtest parameters involved in this simulation are shown in 
Figure [6} 



Parameter 


Value 


K 


1.0 


e 


2 


a 


0.01 


Nmc 


10000 



Figure 6: Backtest parameters 
Let us denote by i9 a and d b the Euler scheme simulation of the compound poisson 



processes i9 a and d b , with dynamics (2.3). Therefore, for a G {a*,a benchmark }, we were 



able to compute the Euler scheme simulation X a (resp. Y a ) of X a (resp. Y a ), starting at 
at t = 0, by replacing r d a (resp. d b ) by d a (resp. ?!r) in equation (2.4) (resp. (2.5)). 

We performed Nmc simulation of the above processes. For each simulation u G 
[1...Nmc] an d for a G {a* , a benchmark } , we stored the following quantities: the terminal 
wealth after terminal liquidation Vjf(u) := L(X a (co), Y a (uj), P(uj)), called "performance" 
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in what follows ; the total executed volume Q total ' a (u)) := X40T] l^t (^0 ~~ ^t-C^)! i an d 
the volume executed at market Q market > a (cu) := ^[ot] ICn(^) a |- Finally, we denote by m(.) 
the empirical mean, by £(.) the empirical standard deviation, by skew(.) the empirical 
skewness, and by kurt(.) the empirical kurtosis, taken over to G [1...Nmc]- 



Quantity 


Definition 


a* 


^benchmark 


Info ratio over T 


m(V T )/a(Vf) 


0.238 


0.104 


Profit per trade 


m{V T )/m(Q total >-) 


1.37 


3.86 


Risk per trade 


a{Vf)/m{Q total '-) 


5.73 


37.21 


Mean performance 


m(Vrp) 


376.08 


773.15 


Standard deviation of perf 




1574.97 


7462.96 


Skew of perf 


skew(Vf) 


0.027 


-0.0468 


Kurtosis of perf 


kurt(Ff) 


3.21 


7.48 


Mean total executed volume 


m(Q total >-) 


274.77 


200.27 


Mean at market volume 


m{Q market >-) 


101.75 


0. 


Ratio market over total exec 


m(Q market '-)/m(Q total '-) 


0.37 


0. 



Figure 7: Synthetis table for backtest. Categories are, from top to bottom: relative per- 
formance metrics, period-adjusted performance metrics, absolute performance metrics and 
absolute activity metrics. 

Figure [7] displayed a synthesis of descriptive statistics for this backtest. We first notice 
that the information ratio over T of a* is more than twice that of a benchmark , Second, 

the per trade metrics can be compared to the half-spread — = 6.25 EUR/contract, and 
we see that although the mean profit per trade is smaller for the optimal strategy, the 
risk associated to each trade is dramatically reduced compared to the benchmark. This 



1400 



1200 



s? 1000 



800 



% 600 

S 
w 

400 



200 



[ 



Opt imal 



-i_ Benchmark 



-15 000 -10000 -5000 5000 

Performance 



10000 



15 000 



Figure 8: Empirical distribution of performance Vf. The graph shows the number of 
occurences for each bin on Nmc = 10000 simulations. 



is confirmed by the empirical distribution of performance, also shown in Figure [8l where 
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the dark blue represents the performance distribution of the optimal strategy and the light 
purple represents the performance distribution of the benchmark strategy. We see that not 
only benchmark has higher standard deviation, but also higher excess kurtosis and heavy 
tails: this is due to the fact that inventory can be very large for the benchmark strategy, 
and therefore it bears a non- negligible market risk (or inventory risk). Finally, we see that 
about 37% of the trades are done with market orders, which indicates that this feature of 
the strategy is relevant when exposed to adverse selection risk (the risk that the mid-price 
moves unfavourably after a limit order execution). 



7 


a(V T ) 


m(V T ) 


6.67 


10- 


04 


13.36 


0.09 


4.44 


10- 


04 


351.16 


20.98 


2.96 


10- 


04 


495.43 


30.97 


1.98 


10- 


04 


649.11 


39.28 


1.32 


10- 


04 


849.05 


106.14 


8.78 


10- 


05 


1048.73 


177.27 


5.85 


10- 


05 


1264.10 


253.50 


3.90 


10- 


05 


1428.10 


309.95 


2.60 


10- 


05 


1546.30 


351.86 


1.73 


10- 


05 


1635.61 


368.15 


1.16 


10- 


05 


1639.65 


332.51 



Figure 9: Varying risk aversion parameter 7: data. 
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Figure 10: Varying risk aversion parameter 7: plot. 



Our last numerical test is devoted to displaying the influence of the risk aversion para- 
meter 7. All other parameters remaining the same, we tested several values of 7 (as 
indicated in Figure [9]) , and characterized the performance of the corresponding strategy by 
the pair (a(Vf), m(V^)), which gives the efficient frontier plot displayed in Figure 10 As 
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expected, a reduction of 7 increases the standard deviation of the strategy: this is due to 
the fact that a small 7 allows for large open position i.e. large inventory, and therefore the 
market risk is greater. For small 7, performance is also better since the investor can sustain 
large inventories, and therefore is less impatient to get rid of it: in particular, the proportion 
of volume executed at market is increasing in 7. We see that in real trading conditions, 7 
must be chosen as small as possible, i.e. the value allowing the greatest performance, but 
maintaining the market risk sustainable. 



5 Best execution problem and overtrading risk 

In this section, we apply our market model framework to a best execution problem. The 
trading objective of the investor is to liquidate Yq > assets over the finite time interval 
[0,T]. She is not allowed to purchase stock during the liquidation period, and may only 
buy back the asset in case of short position. In this context, the investor posts continuously 
a limit sell order (with a volume much larger that the required quantity Yq) at the best 
ask price, and also runs market (sell) orders strategy until she reaches either a negative 
inventory or the terminal date. By doing so, she hopes to trade as much as possible at the 
ask price, and therefore avoiding to cross the spread. 

Mathematically, this means that the investor uses a subset Ae of strategies a = {a make = 
(L a , L b ),a take ) in A such that: 



f (1,0) fort<r, 
j (0, 0) for t > t 

a take = (r„, Cn)n U (r, -Y T ), with r n < r, ( n < 0, 



where r = inf{£ > : Yt < 0} A T. The value function associated to this liquidation 
problem is then defined by 

r f T 1 

v e (t,x,y,p) = sup E tw L(I T ,7 T ,F T )-7 / Y^g(P s )ds, (5.1) 
aeAt 1 Jt J 



for (t,x,y,p) 6 [0, T] x R x P. With the notation in (3.6), the operator corresponding to 
the limit order in Ai is given by C ' = V + r a , while the impulse operator associated to 
the market order in At is defined by: 

Me ip(t,x,y,p) = sup ip(t,x-ep-\e\(- + e) -e ,y + e,p) 

ee [-Ij/1,0] 1 



The dynamic programming equation associated to (5.1) takes the form: 

min [ - - Vv e - T a Vi + 75 , v t - M e v e ] = 0, on [0, Tjxix (0, 00) x P, 
together with the terminal and boundary conditions: 

vi = L, on ({T} xRxtxP)u([0,T) xlxL x P) . 
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The above boundary condition for nonpositive inventory is related to the overtrading risk, 
which is the risk that the investor sold too much assets via the (oversized) limit order at 
the best ask price. This risk occurs typically in execution problems on pro-rata limit order 
book, see jB]. 



Again, in the Levy case (3.11), the value function V£ is reduced into: 



v e (t,x,y,p) = L(x,y,p) + w e (t,y), 
where W£ is solution to the integro-variational inequality: 

{ 

POD 

-A a / 

JO L 2 



Jo 

wi{t,y) 



dwe 2 
nun <; - — -yc p 

$ s 

m(t, y-z)~ W(t, y) + z- + (- + e)(|y| -\y- z\) fi a (dz) ; 



sup w e (t,y + e) - (- +s)(\y + e\ + \e\ - \y\) - e ) 

ee[-y,0] i} 



for (t,y) £ [0, T) x (0,oo), together with the terminal and boundary conditions: 
w e (t,y) = 0, V(t, !/ )e({T}xR)u([0 ) r)xl_). 
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Figure 11: Numerical results for the simple liquidation problem (for cp = 0). On the left 
side, level lines are indicated for the value function vj£. On the right side, numbers indicated 
on the figure represent the quantity to sell in the optimal market order control. 



The associated numerical scheme reads now as follows 

we(tN,y) = 0, 
We(t k ,y) = 0, 

Wi(t k ,y) 



max 



y Ef, 

fc = 0,...,iV-l ) y <0, 

T' 
'I 



(t,y,<p) ; M e ' Y ' {t,y,tp) 



k = 0,...,N-l,y€Y+, 
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where Y+ = Y M n K+, 

7f AY,M (t,y,(p) = <f(t,y) - h^gy 2 + hyc P 

+ X a h( [°° [(p(t,PToj M {y-z))- V (t,y)]fi a {dz) 



o 



+ 



1°° [ 6 ^ + ( 6 - + e)(\y\-\y-z\)]^(dz)) 



and 



M h t A ^ M {t,yM 

5 

= sup [ip(t,Proi M (y + e)) - (- +e)(\y + e\ + \e\ - \y\) -e ]. 

eeY M n[-j/,0] z 

In this case, the optimal policy shown in Figure [TT] is simple to describe. The state space 
is delimited in two zones: when the inventory is small, the HFT must wait for her limit sell 
order to be executed; and when the inventory is large, the HFT must send a market sell 
order to avoid the market risk related to holding a large position. 

The frontier between the two zones (indicated in bold red in Figure [TT| ) can be inter- 
preted as an optimal trading curve, a concept that is extensively documented (see e.g. 0) 
in the optimal execution literature. The optimal trading curve is the inventory that the 
investor should hold, seen as a function of time, in order to minimize overall trading costs. 
Therefore, in the typical setting, the execution strategy consists in trading via market or- 
ders to get as close as possible to the optimal trading curve. Similarly, in our case, we can 



see on Figure 11 that the optimal strategy will behave similarly for large inventories (i.e. 
when above the trading curve): indeed, we observe that the quantities to sell are such that 
the market orders strategy would keep the inventory close to the optimal trading curve, 
if no limit orders were allowed. Now, in our case, we observe two specific features of the 
optimal strategy: 1) the optimal trading curve does not reach at maturity, and therefore 
the HFT has to get rid of her inventory at market at final date to match the constraint 
Yt = 0. This is due to the fact that a supplemental gain is always achievable when the 
limit order is executed. Therefore, this features leads to an execution strategy where the 
final trade is bigger than intermediary trades; 2) below the optimal trading curve, i.e. in 
the region where the HFT trades via limit orders only, the sell limit order is always active, 
and can lead to an execution. Therefore, the inventory is always below the optimal trad- 
ing curve, and the distance between the current inventory and the optimal trading curve 
equals the volume executed via limit orders. This differs from classic pattern-based best 
execution strategies, for example the U-shaped execution strategy that consists in trading 
a large quantity at the beginning and at the end of the liquidation, and trade regularly 
small quantities in between. Indeed, the optimal strategy does not provide a fixed pattern 
for every execution, but provide the optimal action to take given the observation of the 
inventory that is still to be sold and the market characteristics as e.g. the mean traded 
volume at ask per second \ a j2 a , or trades volume distributions at ask /i a . 

Finally, let us notice that this strategy can be interpreted as a convenient way to avoid 
the cost of crossing the spread during the liquidation of a portfolio, but we did not take into 
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account the impact of the market order on the transaction price. In the case of a pro-rata 
microstructure, available volumes offered at best prices are usually about 200 times larger 
than the mean volume of market orders (see [B]), and therefore it is consistent to consider 
that there is no impact on the price for our market orders. Yet, the model can easily be 
modified by adding an impact component in the obstacle operator Mg to take care of this 
effect. We also did not model the possibility that the intensities A a and X b of execution 
processes may vary, and postpone this investigation for future research. 
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